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Summary. — We discuss the high-temperature superconductivity in copper oxide 
ceramics. We propose an effective Hamiltonian to describe the dynamics of elec- 
trons or holes injected into the copper oxide layers. We show that our approach 
is able to account for both the pseudogap and the superconductivity gap. For the 
hole-doped cuprates we discuss in details the underdoped, optimal doped, and over- 
doped regions of the phase diagram. In the underdoped region we determine the 
doping dependence of the upper critical magnetic field, the vortex region, and the 
discrete states bounded to the core of isolated vortices. We explain the origin of the 
Fermi arcs and Fermi pockets. Moreover, we discuss the recently reported peculiar 
dependence of the specific heat on the applied magnetic field. We determine the 
critical doping where the pseudogap vanishes. We find that in the overdoped region 
the superconducting transition is described by the conventional d-wave BCS theory. 
We discuss the optimal doping region and the crossover between the underdoped 
region and the overdoped region. We also discuss briefly the electron-doped cuprate 
superconductors. 

PACS 74.20. -z, 74.72.-h, 74.72.Gh, 74.72.Ek-. 



1. Introduction 

Since the discovery of superconductivity in the copper oxide ceramics (cuprates) [1] , 
both theoretical and experimental efforts have been made to investigate the physical 
mechanism responsible for superconductivity with an unprecedented high transition tem- 
perature T c (for recent reviews see Refs. [2, 3, 4, 5, 6, 7, 8, 9]). Despite this, an under- 
standing of the superconducting pairing of the high-temperature superconductors is still 
lacking. 

Very soon after the discovery of the high temperature superconductors, it was realized [10] 
that superconductivity were intimately related to the square planar CuOi lattice. More- 
over, the physics of the CuC>2 planes were well described by the nearly half-filled Hubbard 
model with moderately on-site Coulomb repulsion. 
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Fig. 1. - Schematic phase diagram of hole-doped cuprate high-temperature superconductors. 
Carrier concentrations 5 near to the maximum critical temperature T c (6 op t) locate the optimal 
doped region. Regions with lower and larger hole concentration are called underdoped and 
overdoped, respectively. The temperature T* is the pseudogap temperature. 



The microscopic model for the description of electrons in the CuOi layers is believed to 
be the effective single-band Hubbard model: ( 1 ) 



(1.1) 



H 



E 

<i.j>y 



4,a Cj 



where c\ a and £j )Cr are creation and annihilation operators for electrons with spin a, U 
is the onsite Coulomb repulsion for electrons of opposite spin at the same atomic orbital, 
and t is the hopping parameter. As is well known [12], the superexchange mechanism 

.2 

yields a Heisenberg antiferromagnctic exchange interaction with J = ^j- between spins 
on the Cu atoms. Since U 3> t, at half-filling the onsite Coulomb repulsion gives an 
insulator where the electron spins are antiferromagnetically ordered. 
The phase diagram resulting from tuning the doping of CuOi planes with holes is shown 
schematically in Fig. 1. At half- filling, corresponding to 6 = 0, we have a Mott insula- 
tor with long range antiferromagnetic order which disappears for temperatures above the 
critical Neel temperature Tn (not shown in Fig. 1). By adding holes into the C11O2 planes 
(6 > 0) the long range antiferromagnetic order is rapidly destroyed. Superconductivity 
occurs beyond a minimal hole content 5 m in — 0.05 — 0.06 where there is no long-range 
antiferromagnetic order. Nevertheless, signs of magnetism persist [13] in the hole-doped 
materials with 8 > 5 m in as local order on a short length scale £af ^ 10~ 6 cm, as observed 
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in muon spin relaxation [14] and neutron scattering [15, 16, 17] measurements. 
The superconducting critical temperature T c rises with increasing hole doping through 
the underdopcd region (5 < S opt ), reaches a peak at optimal doping 6 ~ 6 opt , falls away 
at higher doping 6 > S opt (overdoped region), and then vanishes at the maximal dop- 
ing 8 ma x — 0.3. This peculiar doping dependence of the critical temperature T c is a 
generic feature of the hole-doped copper oxides. Another generic feature of the high- 
temperature cuprate superconductors is the presence of a pseudogap ksT* (ks being 
the Boltzmann's constant) in the underdopcd region S < S op t- The pseudogap was first 
detected in the temperature dependence of the spin-lattice relaxation, Knight shift, and 
magnetic susceptibility studies ( 2 ). The pseudogap temperature T* is much larger than 
the superconducting critical temperature T c in the underdoped region. Moreover, T*(S) 
decreases upon increasing the doping and seems to merge with T c {6) for 5 ~ S op t- In the 
overdoped region, however, the coexistence of the superconducting gap with the pseudo- 
gap is still subject to an intense debate. 

The challenge is to find a theory that is rooted in the microscopic physics of the copper 
oxides which is able to account for the physics of the pseudogap and the supercon- 
ductivity. One of the problem in constructing a satisfactory microscopic theory of the 
high-temperature superconductivity is to isolate that part of the interaction which is 
responsible for both the pseudogap and the superconductivity in cuprates. The purpose 
of this paper is to present an approach which, according to the spirit of the conventional 
superconductivity theory of Bardeen, Cooper, and Schricffer [18] (BCS), allows to iso- 
late an effective Hamiltonian governing the dynamics of the electrons or holes injected 
into the undoped copper-oxide planes. We shall assume that the superconductivity orig- 
inates into the CuOi planes, and thereby we shall neglect the motion along the direction 
perpendicular to the planes. We further assume that electron dynamics in the CuOi 
planes is described by the single-band effective Hubbard model with Hamiltonian given 
by Eq. (1.1). Since the majority of high-tcmpcrature cuprate superconductors are hole- 
doped, we shall focus mainly on hole-doped cuprates. Nevertheless, in Section 6 we 
briefly discuss also the electron-doped cuprates (for a recent review, see Rcf. [19]). 

The plan of the paper is as follows. In Sect. 2 we set up the effective interaction 
between two holes in an antiferromagnetic background and we discuss the main assump- 
tions of the paper. Section 3 is devoted to a general discussion of the physical structure 
of the superconducting ground state. In Sect. 3T we discuss the real-space d-wave bound 
states of two holes; in Sect. 3'2 we address the problem of the superconductivity and 
pseudogap in the underdoped region. In Sect. 3'3 we discuss the d-wave BCS ground 
state and argue that it is relevant in the overdoped region. In Sect. 3'4 we present our 
understanding of the phase diagram of hole-doped cuprates. The physics of the under- 
doped region is presented in Sect. 4 where we discuss the real-space bound states in 
an external magnetic field (Sect. 4T), the upper critical magnetic field (Sect. 4'2), the 
vortex region (Sect. 4'3), the physics of Fermi arcs and quantum oscillations (Sect. 4'4). 
Section 5 is devoted to the discussion of the physics in the overdoped and optimal doped 
regions together with the crossover between the underdopcd region and the overdoped 
region. In Sect. 6 we briefly address the problem of superconductivity and pseudogap in 
the electron-doped cuprates. Finally, in Section 7 we summarize the main achievements 
of the paper and we draw our conclusions. 



( 2 ) See Ref. [8] and references therein. 
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2. The effective Hamiltonian 

In this Section we set up our effective Hamiltonian for the dynamics of holes in the 
nearly half-filled single-band Hubbard model. As discussed in the Introduction, we are 
assuming that the superconductivity in cuprates is related to the square planar CuOi 
lattice. The single-band Hamiltonian is given by Eq. (1.1), where < i, j > denotes a pair 
of nearest-neighbor sites of a square lattice with lattice constant ao (see Fig. 2). Let M 
be the number of sites of the copper lattice and N the number of electrons. In the half 
filled lattice TV = M, i.e. there is exactly one electron per site. In the underfilled case 
N < M there are holes with doping fraction (number of holes per site) : 



(2.1) 



5 = 1- 



N 
M 



In the subspace of the state space of avoided double occupancies of Cu sites, the one-band 
Hubbard model may be canonically transformed into the so-called t — J Hamiltonian: 



(2.2) 



H = —t 



E 

<l,j>,<7 



\& c- 



r T c 



E 

<i,j> 



Si 



* 1 . - 
Sj - - n, u. 



where 



(2.3) 



i,a (1 



J = 



4^ 
U 



and S% is the spin-1 operator at site i. ( 3 ) Since at half-filling hi t - a = 1, 
Hamiltonian Eq. (2.2) reduces to the square lattice Heiscnbcrg Hamiltonian: 



the t — J 



(2.4) 



H 



J ^ Si ■ Sj 



<i,3> 



The ground state of the latter Hamiltonian is known to be the antiferromagnetic Neel 
state. When the system is doped with holes away from half filling (5 > 0), the dynamics 
is no longer described by the Heisenberg Hamiltonian Eq. (2.4), but one must use the 
t — J Hamiltonian Eq. (2.2) with the supplementary constraint of no double occupied 
lattice sites. It is easy to see that the motion of dilute holes in the antiferromagnetic Neel 
background is strongly frustrate. Indeed, when a hole hops from one site to its nearest 
neighbor, a spin is also moved onto one lattice site with the wrong spin orientation 
with respect to the antiferromagnetic background. On the other hand, a pair of holes 
with nearest-neighbor distance (see Fig. 2) can move almost freely without destroying 
the antiferromagnetic background. The effective Hamiltonian for the propagation of the 
holes in the antiferromagnetic background can be written as [21, 22]: 



(2.5) 



H 



U 



E E ^'h^+ ia ° + i a o)^h{r) 



• ■j 



( 3 ) See, for instance, Ref. [20] and references therein. 
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Fig. 2. - The idealized square copper plane with lattice spacing ao. The dashed ellipses represent 
two holes in an antiferromagnetic background. 



where ip h (r), iphir) are creation and annihilation operators for holes at the lattice site 
r and the sum over the direction vectors i and j is restricted to next-nearest neighbor 
lattice sites. Note that in the effective Hamiltonian Ho we do not display the spin of the 
holes since the antiferromagnetic background forces the two holes to have antiparallcl 
spins. In other words, the motion of a hole must not disturb the antiferromagnetic 
background. Moreover, in this antiferromagnetic background approximation, it turns 
out that there is an effective attractive two-body potential between nearest-neighbor 
holes [21]. However, it has been pointed out [23] that the motion of a pair of holes, 
which naively should be mobile, is frustrated due to fermionic nature of the background 
spins. Indeed, for states with two holes as nearest neighbors in an antiferromagnetic 
background it is easy to see [24] that there are configurations that differ by the exchange 
of two electrons. The destructive interference of these configurations leads to holes with 
very large effective mass. This problem can be avoided if we consider pairs of hole with 
distance much greater than the lattice spacing. Accordingly, we shall suppose that two 
holes with distance ao *C vq -C £af, where £af is the antiferromagnetic local order 
length scale, are subject to an effective attractive two-body potential. Our proposal is 
quite similar to the spin-bag theory [25] where the pairing is due to a local reduction of 
the antiferromagnetic order (bag) shared by two holes. 

Thus, we are lead to consider the following reduced interaction Hamiltonian: 



(2.6) 



H, 



= § / dridr 2 W h>t (ri) ^4(^2) V(n - r 2 ) ^h,i(r 2 ) VvKn) • 
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The simplest choice for the two-body potential V{f{ — fz) is: 

oo f[ = r 2 

(2.7) V(fi-f 2 ) = { -V \n-f 2 \ < r (S) 

otherwise 

where the hard-core at the origin is imposed to avoid configurations where two holes are 
on the same site. Moreover, since we are considering the limit of low carriers, S <C 1, the 
hard-core potential does not allow two holes to come too close, thereby implementing 
the no-double occupation constraint. As concern the range of our potential ro(<5), it is 
natural to choice ro(S) of the order of the coherence length £o — 6ao- So that we assume: 

(2.8) r (S) = 6a (l - pj* , S c ~ 0.35. 

The peculiar dependence on the doping fraction 6 arises as follows. The attractive two- 
body potential originates from the interaction of two holes in the antiferromagnetic back- 
ground which extends over a distance of order £af- Since the injected holes tend to 
destruct the antiferromagnetic order, the area of the antiferromagnetic islands should 
decrease with increasing 8. In addition, as discussed in the Introduction (see Fig. 1), the 
superconducting instability disappears at the maximal doping S max ~ 0.30. The simplest 
parameterization which takes care of these effects leads to our Eq. (2.8). Regarding the 
parameter Vq in Eq. (2.7), we have fixed this parameter such that there is at least one 
real space d-wave bound state. To this end it suffices to assume: 

oy-2 

(2.9) V Q ~ 2 J = ^- . 

In the following we will use the numerical values: 

(2.10) a ~ 4.0 10^ 8 cm , t ~ 0.11 ev , U ~ 1.1 ev , 

that are appropriate for a typical cuprate. Note that the resulting Neel temperature: 

(2.11) T N ~ — - lO 3 ^ 

is of the order of the observed Neel temperatures in cuprates. 

Before concluding this Section, we would like to rewrite the Hamiltonian Ho Eq. (2.5) in 
the limit of low-lying excitations. Observing that H can be diagonalized by writing: 

(2.12) 4(r) = -L £ exp (ik ■ r) $ h {k) , 



so that: 



(2.13) i> h {k) = -j= ^expHfc-r) ^ h (f) 
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We get: 

(2.14) H = i>l(k)i> h (k), 



where: 



2t 2 

(2.15) £j = — [cos(2k x a ) + cos (2k y a ) ] 



We now take the small-k limit: 

(2-16) H = M)Mk) 

h 

with: 

( 2 - 17 ) m *h = ■ 

8 TJ a 

So that we may write for the effective Hamiltonian: 

(2-18) H = f&$l<r(f) f^^W(r), 

a=U J \ Zm h J 

where we make explicit the spin degrees of freedom. Using the numerical values of the 
parameters we obtain: 

(2.19) m* h ~ 5.41 m e , 

where m e is the electron mass. To summarize, our effective Hamiltonian for low-lying 
excitations may be described by: 

(2.20) H eff = Ho + H in t , 

with Ho and Hint given by Eqs. (2.18) and (2.6) respectively. It is worthwhile to stress 
that the arguments developed in this Section cannot be considered as a truly microscopic 
derivation of the effective Hamiltonian Eq. (2.20). 

3. The superconducting ground state 

For later convenience it is useful to think about the wave function of the superconduct- 
ing ground state within the BCS theory. Let (p(ri —r^) be the wavefunction of a Cooper 
pair [26] , thus the superconducting ground state may be expressed as an antisymmetrized 
product of identical pair functions: 



(3.1) Q(fi, ...,f N ) = ]T (-1) P PMn - f 2 ) • • • ^(fjv-j - r N )} , 

p 
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in which the sum is over all permutations P of the N particles. The expression Eq. (3.1) 
may be used to point out the analogies and differences between the superconducting 
ground state and that of a Bose-Einstein condensation of pairs [27, 28] . Indeed, a Bose- 
Einstein condensate of pairs would also be given as a product of identical pair wave- 
functions. The process of antisymmetrizing would not change the character of the state 
very much if the size of the pair wavefunction, £o; were small compared with the average 
spacing between pairs. In this case, wavefunctions in the sum differing by the exchange of 
single members of two or more pairs would not overlap very much, and the system would 
behave qualitatively like a condensed Bose-Einstein gas. However, antisymmetrization 
makes a major difference if the size of the pair wavefunction is large compared with 
the spacing between pairs. In this case, it is convenient to rewrite the superconducting 
ground-state wavefunction Eq. (3.1) in the second quantization formalism. To this end, 
we introduce the creation operator of a pair at rest: 



&t = / dfdf <p(f - f 2 ) tp\(f{) ^\{f^) 



(3.2) 

Then, the wavefunction Eq. (3.1) can be written as 



(3.3) 



m > 



N 



N'\ 



7T l°>' 



where |0 > is the vacuum, and N' — £ is the number of pairs. Now, from the identity [29] 

1 r 2 * 

(3.4) \Q >N = — d< p e -' lN * exp [e^tf] |0 > , 

2tt Jo 

we argue that the state: 
(3.5) 



\n>Bcs= exp^&t] |0> 



is the superconducting ground state for a system with varying number of particles. Obvi- 
ously, in that case the relevant Hamiltonian is H — fj,N, where is the chemical potential. 
From Eq. (3.2) we obtain: 



(3.6) 



where: 
(3.7) 



-f 2 ) = y^exp [if (fx - r 2 )} <p{q) 



Using Eq. (3.6) it is easy to find: 

(3.8) |n >bcs = n l 1 + ^ V® 4,t 6 -?,J 1° > > 
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with the normalization factor Af$ . Indeed, Eq. (3.8) agrees with the well-known BCS 
variational ground state: 



(3.9) 



\o >bcs = n [ u ® + v ^ 4,t £ -L 



10 > 



with the normalization: 
(3.10) 



Hq)\ 2 + WW 



1 



As we discuss later on, the superconducting ground-state wave function in the form of 
Eq. (3.1) will be of relevance in the underdoped region, while in the overdoped region 
the BCS variational ground state Eq. (3.9) would be more appropriate. 

3'1. Real space d-wave bound states. - We are interested in the problem of two holes 
in interaction with the effective two-body potential Eq. (2.7). The resulting Schrodinger 
equation is: 



(3.11) 



2ml 



[Vl + V| 2 ) + V(n - f 2 ) 



®{r\ - r 2 ) = E $(f=i - f 2 ) 



According to the reduced interaction Hamiltonian Eq. (2.6) the holes have antiparallel 
spins, so that the wavcfunction <S>(ri — must be symmetric under hole interchange. 
Introducing: 



(3.12) 

we rewrite Eq. (3.11) as: 
(3.13) 



2ml 



2Vi 



ri - r 2 



iv 2 - 

2 R 



R 



V(r) 



n + f 2 



$(R,r) = E <f>(R,r) , 



whereupon: 
(3.14) 



$(R,r) = exp (iK ■ R) <p(r) 



(3.15) 



E 



h 2 k 2 

4mJ 



(3.16) 



(p(r) = — A ip(r) 



Without loss in generality, we may assume that the center of mass of the pair is at rest, 
P = HK = 0. Since the potential is central we may adopt polar coordinates. We have: 



(3.17) 



cxp(im6») 
Vyr,0) = 7= — <^m(r) , 
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Inserting into Eq. (3.16) and using Eq. (2.7) we get: 
(3.18) <^ m (r = 0) = 0, 



(3.19) 



dr 2 



1 d 
r dr 



+ -± (-A m + V ) 



ip m (r) = , < r < r (S) , 



(3.20) 



dr 2 



1 d 

r dr 



h 2 



■Ar, 



Vm{r) = , r (6) < r 



Since (p m (r) is even and must satisfy Eq. (3.18), we have m = 2,4, .... With our choice 
of the parameters, Eqs. (2.8), (2.9) and (2.10), it is straightforward to check that there 
are bound states (A m > 0) only for m = 2 (d-wave). In this case we get: 

(oo-n / r \ _ / AJ 2 (kr) r < r (S) 

(3 ' 21) V2[r) - \ BK 2 (kr) r (6) < r 

where Ji(x) and K2{x) are Bessel functions, 
(3.22) 



m 

H 2 



A 2 , fc 2 = ^ (V - A 2 ) 



and 



(3.23) 



Hr (S) 



K' 2 (kr (S)) 
K 2 (kr (S)) 



kr (6) 



j'2(kro(S)) 
J 2 (kr (S)) 



With the further definitions: 



(3.24) 



xq = k 7- (5) 



m h V o 

h 2 



' c = F' 



one can rewrite Eq. (3.23) as: 



(3.25) 



if 2 (x - C 2 ) 

K 2 (x ^W 2 ) 



£o C 



J 2 { x o 
■72 (xo C) 



It is not too hard to see that Eq. (3.25) allows non trivial solutions if the effective 
wall-depth xo satisfies: 



(3.26) 

where xq is the solution of: 
(3.27) 



xo > xq , xo — 3.8317 



J 2 (x) 



Mx) 
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Fig. 3. - The gap A 2 (<5) versus the hole doping fraction 5. According to Eq. (3.28) S* ~ 0.207. 



From Eq (3.26) it follows that there are non trivial solutions for hole doping such that: 



(3.28) 



< S < 5* 



S* 



0.207 



In fact, in Fig 3 we display A2(<5) versus the hole doping fraction 5. We see that the 
gap A2(S) decreases rapidly with 5 and vanishes at S = S* according to Eq. (3.28). Note 
that, with our parameters, we have A(0) ~ 486 K ■ ks which should be compared with 
Vq ~ 1021 K ■ k,B- From Figure 4, where we display the normalized wavcfunction v?2(?")j 
Eq. (3.21), for a typical value of 5, we see that our bound-state wavcfunction extends up 
to distance r ~ 6ao- 

Since the eigenvalue equations Eqs. (3.19) and (3.20) depend on m 2 , the general d-wave 
wavefunction is a linear combination of tp2(r) and <p-2(r). From the geometry of the 
C11O2 planes we are lead to assume: 



(3.29) 



Mr, 8) = Mr) - cos (20) 



where the coordinate axis are directed along the Cu — O bond directions as in Fig. 2. This 
is the most natural choice since the wavefunction is sizable along the Cu — O bonds and 
vanishes at 

the wavefunction Eq. (3.29). We have: 



±?. For later convenience, we need to evaluate the Fourier transform of 



(3.30) 



<P2(k) = / dfip 2 {r,9) exp(ifc-r) 
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Fig. 4. - The normalized d-wave wavefunction f>2 for S = 0.18. 



r/a„ 



Assuming that the k x , k y axis are oriented along the copper-oxygen bond directions, and 
using the expansion: 



+oo 



(3.31) exp (ik ■ r) = cxp (ikr cos 9k r ) = exp [im(0fc r + — )] J m (fcr) , 



rn— — oo 



we readily obtain: 

(3.32) 

where: 



(p 2 {k,e k ) = <p 2 {k) cos(20 fe ) 



(3.33) 



Vi{k) = -27T / dr r </2 2 (r) J 2 {kr) 
Jo 



Remarkably, we see that the Fourier transform of the wavefunction vanishes along the 
nodal directions: 



(3.34) 



fc = ±T 



while it is sizable along the antinodal directions: 



(3.35) 



0, ± 
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These features are reminiscent of the observed pseudogap in the underdoped region. 
Therefore, we are lead to identify ^2(8) with the pseudogap and to introduce the pseu- 
dogap temperature according to: 

(3.36) k B T*(S) = l -& 2 (8) . 

3'2. Superconductivity in the underdoped region and the pseudogap. - We have shown 
in the previous Section that, as long as 8 < 5*, it is energetically favored to pair two 
holes into a d-wave bound state with wavefunction given by Eq. (3.29). Then, we may 
construct the superconducting ground state according to Eq. (3.1) with ip(r\ — fj) = 
Vi^Ti ~ Fj). However, although the formation of Cooper pairs is essential in forming 
the superconducting state, its remarkable properties, i.e. zero resistance and perfect 
diamagnetism, require phase coherence among the pairs. In the underdoped region <5 <C 
6* we infer from Fig. 4 that the size of the pair is small with respect to the average distance 
between holes. So that the ground state wavefunction Eq. (3.1) reduces to the ground 
state of a Bose-Einstein gas of paired holes. In this case the required phase coherence 
is established by the condensation of pairs. However, in two spatial dimensions it is 
well known that there is no Bose-Einstein condensation. Thus we are lead to conclude 
that in the extremely underdoped region 5 <C S* there is a pseudogap A2(<5) but not 
superconductivity. However, by increasing S it will be a certain minimal value of the 
hole doping fraction 5 m i n where the pairs become to overlap. To estimate 8 m i n , we 
observe that the pair wavefunction <p2(f) extends up to distance r ~ £0 ~ 6oq. Thus, to 
have overlap between different hole pairs we need at least two holes in a square of size 
6ag. Thereby we estimate: 



(3.37) S min ~ p ~ 0.056 



which agrees, remarkably, with the observed value S m i n ~ 0.05 — 0.06 for the onset of 
the superconductivity in the hole-doped cuprates. Once the pairs begin to overlap, the 
phases of the pairs are locked to a constant value. Indeed, it is well known ( 4 ) that a 
condensate with varying phase Q(R) contribute to the ground state energy with: 

(3.38) H e = ^ J dR\Ve(R)\ 2 , 

where K s is the so-called phase stiffness. In our case we have: 

h 2 



(3.39) K s 
with the supcrfluid density given by: 



2ml 



(3.40) n s ~ — — , S min < S . 

Za 

As is well known [31, 32], the phase coherence of the condensate survives up to the 



( 4 ) For instance, see Ref. [30]. 
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Fig. 5. - The pseudogap and Berezinskii-Kosterlitz-Thouless temperatures (in Kelvin) versus S. 
We assume 8 m in ~ 0.05. 



Berezinskii-Kosterlitz-Thouless (B-K-T) critical temperature: 

(3.41) k B Tb-k-t | K s {T B -k-t) . 

For temperatures above Tb-k-t the phase coherence and the superconductivity is lost 
due to the thermal activation of vortex excitations. If we neglect the temperature de- 
pendence of the phase stiffness, we obtain: 

(3.42) k B Tb-k-t — % — r -j — T — 6 , S min < 8 . 

8 m\ al U 

Using our numerical values for t and U we get: 

(3.43) Tb-k-t - 401 K 8 , S mm < S . 

In Figure 5 we display the Bcrezinskii-Kostcrlitz-Thouless critical temperature together 
with the pseudogap temperature T* versus the hole doping fraction by assuming 5 m i n ~ 
0.05. We see that in our approach the superconductivity transition in the undcrdoped 
region is described by the B — K — T transition at T c = Tb-k-t, which coexists with the 
much higher pseudogap temperature T*. This picture is supported by a large amounts 
of observational data which, however, will be discussed in details in Sect. 4. 

3'3. d-wave BCS ground state. - To complete the phase diagram of the hole-doped 
cuprate superconductors we must address the problem of the ground state for 8 > S*. 
Since the pseudogap A2(<5) vanishes at 8 = 8*, we do not have real-space paired holes 
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and, whence, the B — K — T critical temperature cannot be defined. On the other hand, 
as long as S < 5 C ~ 0.35, the attractive two-body potential Eq. (2.7) is short-range and 
can be considered as a small perturbation. As a consequence, in the overdopcd region 
S > 6* the normal state is described by the Hamiltonian Ho, Eq. (2.18). In other words, 
the normal state properties are those of ordinary Fermi liquid metals that is characterized 
by hole quasiparticlcs with effective mass given by Eq. (2.17) and density: 

(3.44) n h ~ ?-t- . 



In fact, in the overdopcd side of the superconducting phase diagram, experiments strongly 
suggest that above the critical temperature the full Fermi surface is restored with well- 
defined quasiparticles. Thereby this confirms the applicability of a Fermi liquid picture. 
Moreover, quantum oscillation experiments [33] in an overdopcd high-temperature su- 
perconductor are consistent with hole quasiparticles with density given by Eq. (3.44) and 
effective mass: 

(3.45) m* h = (4.1±1.0)m e 
in fair agreement with our estimate, Eq. (2.19). 

The superconductive instability is driven by the short-range attractive interaction be- 
tween the quasiparticles. We see, then, that the pairing is in momentum space, so that 
the relevant superconducting ground state is the BCS variational ground state Eq. (3.9) 
with Hamiltonian H — fiN, where the chemical potential is essentially the Fermi energy: 

(3.46) (J, ~ e F , e F = — -f , k F ~ y/2Tm h . 

The relevant gap equation has been discussed since long time [34] : 

-» 1 V(k-k') A(fc') 

(3-47) A(k) =--E 



k' 



q, + |A(fc' 



where: 



(3-48) ^ = ^— - e F 



h 2 P 



and 

(3.49) V(k - k') = I dr cxp [-i(k - k') ■ r\ V(r) 



Using the expansion Eq. (3.31) we may recast Eq. (3.47) into: 

- +oo 

(3.50) A(k) = exp(im^ 



V m (k,k')A(k') 



m— — oo 



e, + iA(fcoi 2 
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with: 

(3.51) V m {k,k') = 2tt 
Equation (3.50) suggests that: 

(3.52) A(k) = A(M fc ) = E^W^W' 



In the weak coupling limit one expects that the mixing of different partial wave £ will be 
small [34]. Therefore Eq. (3.50) reduces to: 

(3.53) A m (k,9 k ) ~ -- 2^ exp(jm6» fc fe^ 



^ + | Am(A ., )0fc)) |2 



Since the BCS gap in the limit of localized pairs is proportional to the bound-state wave 
function in momentum space, we conclude that Eq. (3.53) admits non trivial solutions 
for the d-wave gap (to = ±2). Moreover, from Eq. (3.32) we infer that: 

(3.54) A 2 (Mfc) = A BCS (k,6 k ) = A BCS (fc) cos(2# fc ) . 

Thus, we obtain the following gap equation: 

V 2 (k,k')A BC s(k') 

^ [CUB \ AU k > )\ , 



(3.55) A BCS (k) * -^[cos(2^)] 



q, + [A B cs(k')cos (26V)] 2 



Since the gap is sizable on the Fermi surface, we may further simplify Eq. (3.55) as: 

dk' [cos(26V)] 2 



(3.56) 



1 



-V 2 



( 2?r ) 2 J & + [A BCS cos(2^)] 2 



where A B cs = Abcs^f) and V<i = l^ikp, kp). A standard calculation gives: 



(3.57) 



h 2 



JL 



2n 



d6 



[cos (26)f 



v/^ + [A scs cos(20)] 2 



where e c is an energy cut-off much smaller than the Fermi energy. Performing the 
integrals and using the approximation: 



(3.58) 
we get: 
(3.59) 



archsinh 



1BCS 



cos (26>)| 



In 



2e r 



Ascsl cos (26>)| 



a 4£c r 1 i 

A BC s - —;= cxp - — 

\/e A 2 
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where: 

m*Vn r r °( s ) 

(3.60) A 2 = -Jp y drr[J 2 (fc^r)] 2 . 
To obtain the critical temperature, we note that [35]: 

(3.61) ^ * , 

7 = 0.577216... being the Euler's constant. As concern the high-energy cut-off, it is 
natural to identify s c with A 2 ((5 = 0) which, indeed, is much smaller than the Fermi 
energy in the range of hole doping fraction of interest. In summary, we obtain for the 
critical temperature T c = Tecs'- 

(3.62) k B Tbcs — A 2 (0) exp[-— ] . 

7T A 2 

Equation (3.62) shows that the BCS critical temperature depends on the hole doping 
fraction. In fact, in Fig. 6 we compare Tbcs(S) with the pseudogap temperature T*(5), 
It is evident that for 5* < S < S max ~ 0.30 we have a non-zero BCS d-wave gap. Note 
that the undulations on the critical temperature curve is an artifact of the sharp cut-off 
of the effective interaction potential at r — tq (S) . 
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Fig. 7. - The pseudogap, Berezinskii-Kosterlitz-Thouless, and BCS critical temperatures (in 
Kelvin) versus the hole doping fraction S. The region enclosed by the ellipse where the three 
temperatures are comparable is the optimal doped region. 



3'4. Phase diagram of hole-doped cuprates. - We are, now, in position of discuss the 
phase diagram of the hole-doped cuprates. In Fig. 7 we display the Berezinskii-Kosterlitz- 
Thouless (Tb-k-t), the pseudogap (T*), and the BCS {Tbcs) temperatures versus the 
hole doping fraction 5. In our approach the superconducting dome-shaped region of 
the phase diagram of oxide cuprates is determined by the region enclosed by the two 
critical temperatures Tb-k-t and Tbcs- The Bcrczinskii-Kostcrlitz-Thouless temper- 
ature increases with the doping 5, while the BCS critical temperature decreases with 
increasing hole doping. These temperatures are comparable for S ~ S op t- Accordingly, 
we individuate three regions in the phase diagram: 



• the underdoped region S m i n < 5 < S opt ; 



• the optimal doped region S ~ 5 op t ; 



• the overdoped region 6 opt < ^ < S max . 

In Fig. 7 we have also displayed the pseudogap temperature T* (S) which decreases rapidly 
with increasing doping, meets the superconducting dome in the optimal doped region, 
and vanishes beyond a critical doping level 6* within the optimal doped region. In the 
following Sections we shall discuss the physical properties which characterize the three 
regions and compare with selected experimental observations supporting our proposal. 
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4. The physics of the underdoped region 

The underdoped region in all hole-doped high-temperature superconductors is char- 
acterized by the presence of the pseudogap ( 5 ). In general it is assumed that the pseu- 
dogap has a d-wave structure. Moreover, there is a general agreement on the doping 
dependence of the pseudogap. For instance, recently the authors of Rcf. [38] used angle 
resolved photoemission spectroscopy to probe the electronic excitations in underdoped 
Bi2Sr2CaCu20s+s- These authors find evidence of a nodal liquid whose excitation gap 
vanishes only at points in momentum space which are consistent with the gap structure 
of a d-wave superconductors. In particular, these authors show that the measured gap 
in momentum space is consistent with: 

(4.1) A 2 (5,e k ) = A 2 (6) cos(20 fe ) , 

where A2(<5) monotonically increase with undcrdoping. In fact, Eq. (4.1) is consistent 
with our Eq. (3.32) and the doping dependence of the pseudogap (see Fig. 3). However, 
it should be stressed that many experiments report a gap that deviates from the simple 
relation in Eq. (4.1) [39, 40, 41, 42, 43, 44, 45, 46]. In fact, several experimental studies 
are consistent with the presence of two gaps, namely the pseudogap and a peculiar nodal 
gap. We believe that this nodal gap is related to the dynamics of the nodal Fermi liquid 
discussed in Sect. 4'4. Then, we see that the nodal gap has naturally a d-wave symmetry 
with gapless nodes. On the other hand, the pseudogap being defined as the energy needed 
to break the real-space two-hole bound state, docs not need to depend on the geometry 
of the Fermi surface. However, a full account on this subject lies beyond the aim of the 
present paper. We plan to present a detailed discussion of this point in a separate paper. 
As concern the pseudogap temperature, in Fig. 8, where we plot jr- versus fc 2 ^?„ for 
various cuprates [47], we compare the experimental data with our relation Eq. (3.36) 
which can be rewritten as: 

(4.2) = 4.0 . 

y ' k B T* 

Indeed, we see that the data are consistent with Eq. (4.2). We also plot in Fig. 8 the 
d-wave BCS relation ( 6 ): 

(4.3) ~ 4.3 . 

v ' k B T* 

Since Eq. (4.3) is quite close to our Eq. (4.2), the data are consistent with both relations. 
In our opinion, this unfortunate circumstance is the origin of confusion and several mis- 
interpretations in the literature. 

In conclusion, according to our previous discussion, the observed gap opening in the 
density of states (the pseudogap) marks the onset of the hole pairing at rather high tem- 
peratures T <T*, but the critical temperature is set by the onset of phase coherence at 
lower temperatures T < T c = Tb—k—t- The above point of view is supported by several 
experiments. We believe that the most convincing evidence comes from Ref. [48]. In 



( 5 ) For a review on the pseudogap in cuprates, see Refs. [36, 37]. 

( 6 ) The d-wave BCS critical temperature will be discussed in Sect. 5. 
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Fig. 8. — ^r- versus £ j for various cuprates compared to our relation ^ g ^, = 4 (dashed line) 
and the d-wave BCS relation fc 2Z ^, = 4.3 (dotted line). The data have been extracted from 
Fig. 5 of Ref. [47]. 



fact, by using ultrathin films of Ca-substituted Y Ba^Cu^O^-^, the authors of Ref. [48] 
found a clear evidence of the Berczinskii-Kosterlitz-Thoulcss transition. In particular, 
they observe a remarkable scaling of the critical temperature with n s in accordance with 
Eqs. (3.39) and (3.41). Further support on this picture comes from measurements of the 
high-frequency conductivity reported in Ref. [49]. 

To summarize, we may conclude that the key features of the underdoped side of the 
phase diagram are controlled by very strong hole pairing that is phase-disordered by 
fluctuations [50]. In the following subsections, from a wealth of possible material, we 
have chosen to discuss some selected topics that we think well illustrate the physics of 
the underdoped region of hole-doped high-temperature superconductors. 

4'1. Real space bound state in magnetic fields. - We are interested in the dynamics 
of hole pairs in presence of an external magnetic field perpendicular to the copper-oxide 
plane: 

(4.4) H= VxA , 1 = (A 1 (r),A 2 (r),0) . 

In this case the Schrodinger equation Eq. (3.11) becomes: 




(4.5) +V(r 1 -r 2 )$(r 1 ,r 2 ) = E<f>(r u r 2 ) 
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(4.6) 



A = --fx H 
2 



where the angular momentum is conserved. Using Eq. (3.12) we recast Eq. (4.5) into: 



4m 
(4.7) 

Thus, we can write: 

(4.8) 

and 

(4.9) 

such that: 
(4.10) 



1 -ihV a - —A(R)\ + -\ (-ihVp - £ A(f S 



2c 



$(i?,f) 



V(r) $(i?,r) = E<f>(R,r) 



$(i?,r) = *(i?)<p(r) 



E — E con d — A , 



4m 



h L 



2e - 
■ifiVa : A(i?) 



*(i?) = £ cond , 



(4.11) 



ml \ 2c 



V(f) 



tp(r) = — A y(f) 



The wavefunction ^(R) describes the condensate. In fact, ^(R) is obtained from $(R, r) 
after averaging over the internal degree of freedom. Moreover, if we write: 



(4.12) 



*(i?) = \^{R)\ cxp[ l e( J R)] , 



we see that in presence of an external magnetic field the condensate acquires a non- 
uniform phase. For an almost uniform condensate we have: 

(4.13) I = constant , 

and Q(R) is a slowly varying function. In this case we obtain: 

(4.14) E cond ~ X - (2ml) v?{R) , 

where the supercurrent velocity is given by the well-known relation [51]: 



(4.15) 



2m* h v s (R) = HV A Q(R) - ^ A{R) 
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Let us, now, consider Eq. (4.11) in polar coordinates: 
(4.16) 



-h 2 f (P_ ld_ 

dr 2 r dr 



1 d 2 



e 2 H 2 
\&m* h c 2 



iehH d 



2m*c 39 



V(r) 



(p(r, 6) = -Aip(r, 



Again, after writing: 
(4.17) 



<p(r,0) 



exp (irrib 



we have: 
(4.18) 
-h 2 ( d 2 Id 



ml \ dr 2 r dr r 2 



e 2 H 2 
16m^c 2 



V(r) 



, . . . , TT . ehmH 
Vm\r) = -(A m (H) + — — — ) tp m {r) 
2m h c 



This last equation shows that, as expected, the external magnetic field lifts the degeneracy 
m — > —m. 



4'2. The upper critical magnetic field H c2 . - We have already seen that for H = 
two holes give rise to a d-wave bound state. In presence of an external magnetic field the 
d-wave bound states are the solutions of our pervious Eq. (4.18) specialized to m = ±2: 
(4.19) 



_ l£ - 1 

m* h V dr 2 r dr r 2 



e'H 2 
ia>m* h c 2 



V(r) 



P±a(r) = -(A ±2 (H) ± ^- r ) ip± 2 (r) . 



We may, now, introduce the pair-breaking critical filed if** such that: 



(4.20) 



A± 2 (h£,6) = 



To estimate the pair-breaking critical magnetic field we note that: 



(4.21) 



A ±2 (0,<5) = A 2 (S) 



where the pseudogap A 2 (S), already discussed in Sect. 3T, can be approximated quite 
accurately by: 



(4.22) 



A 2 (S) ~ A 2 (0) 



(-) 1 ' 5 



A 2 (0) ~ 41.91 mev 



Moreover, it is easy to see that the Zeeman term can be neglected. In fact, using 
Eq. (2.19) wc find: 



(4.23) 



ehH 



mTc 



2.14 1CT 5 ev H(T) , 
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where H(T) means that the strength of the magnetic field is measured in Tesla ( 7 ). 
Accordingly, we get from Eq. (4.19): 

(4.24) A ±2 (H,S) * A 2 (S) - < r 2 > . 

Since < r 2 >~ £ 2 , we obtain: 

(4-25) H£ * ^^A 2 (J) , 

which leads to: 



(4.26) 




Equation (4.26) shows that the pair-breaking critical field is very high. This means that 
the pseudogap is not affected by applied magnetic fields H < 10 2 T in accordance with 
several experimental observations. 

In conventional superconductors the pair-breaking critical filed is of the same order of 



( 7 ) IT = 10 4 G. Even though we are using CGS units, it is customary in the literature to 
express the applied magnetic field in Tesla. 
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the depairing critical field. In our case, it turns out that the depairing critical field is 
much smaller than the pair-breaking critical field. As a consequence, the upper critical 
field H C2 is determined by the depairing field. It is worthwhile to recall that in type-II 
superconductors in the mixed state (Schubnikov phase) between the lower critical field 
H Cl and the upper critical field H C2 , the magnetic flux penetrates in a regular array of 
flux tubes, each carrying a quantum flux: 

(4.27) <fo = =^p ~ 2.07 1(T 7 Gem 2 . 

Within each unit cell of the array there is a vortex of supercurrcnt concentrating the 
flux toward the vortex center. To estimate the critical velocity v c we use Eqs. (4.9) and 
(4.14) to get: 

(4.28) E cond ~ m* h v 2 ~ A 2 (H,5) ~ A 2 (<5) . 
Accordingly, the critical supercurrent is: 



(4.29) J c = 2 e n s v c ~ e — v c ■ 

From the Maxwell equation V x H = — J we easily get the following estimate 



Since the magnetic field can be considered almost uniform over distance of order of the 
London penetration length A, and considering that there are N vort ~ (^) 2 = k 2 vortices 
in a region of area A 2 , we obtain our final estimate of the upper critical field: 



(4.31) H C2 {5) 



2 4tt eS A 2 (S) 



In Fig. 9 we display our estimate of the upper critical field Eq. (4.31) in the underdoped 
region 6 < 5*. Wc see that in the deep underdoped region the upper critical field 
decreases with underdoping in accordance with recent measurements [52, 53, 54]. In 
Fig. 9 we also display the depairing critical field as determined from the Nernst signal 
in Refs. [52, 53]. It is known that the cuprate superconductors have a large Ginzburg- 
Landau parameter k ~ 10 2 and are, therefore, extreme type-II superconductors. In fact, 
to compare quantitatively Eq. (4.31) with experimental data we choose k ~ 130. It is 
gratifying to see that the doping dependence of our critical magnetic field as implied by 
Eq. (4.31) seems to follows quite closely the experimental data. 

4'3. The vortex region. - In the underdoped region of the high-temperature cuprate 
superconductors there is convincing evidence for vortices at temperatures significantly 
above the critical temperature [52, 53, 55, 56]. Indeed, our previous determination of 
the upper critical magnetic field shows that H C2 (5) is different from zero even for 5 < 
5 m i n ~ 0.05, i.e. outside the superconducting dome in the phase diagram. To strengthen 
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Fig. 10. - The pseudogap, Berezinskii-Kosterlitz-Thouless, and vortex temperatures (in Kelvin) 
versus S. The vortex region lies in between Tv and Tb-k-t- 



this point, following Ref. [57], we define a critical vortex temperature Ty{8) by the 
phcnomcnological relation: 

(4.32) k b T v {8) ~ 2.1 /is H C2 {5) 

where \xb is the Bohr magneton. The meaning of Ty{8) is that vortices are absent for 
T > Ty{5). In Fig. 10 we compare the vortex temperature to the critical Berezinskii- 
Kosterlitz-Thouless and pseudogap temperatures. For temperatures Tb~k-t{5) < T < 
Tv(S) (the so-called Nernst region) one should observe a large Nernst signal well above 
the superconducting critical temperature, in accordance with observations. 
We would like to conclude this Section by addressing the puzzling observations of dis- 
crete core bound states in isolated vortices. In a classical paper [58], it was predicted 
the existence of electronic states bound to the core of an isolated vortex at energies 
below the superconducting gap in the case of s-wave BCS superconductors. Since the 
spacing of the discrete spectrum turns out to be very small, these bound states lead to a 
continuum which manifests itself as a zero-bias conductance peak in scanning tunneling 
spectroscopy measurements. In high-temperature superconductors the scanning tunnel- 
ing spectroscopy of isolated vortices revealed striking electronic signatures which cannot 
be understood within the conventional BCS s-wave or d-wave theory ( 8 ). Indeed, the 
vortex cores give evidence of a pair of states with energy which scales with the supercon- 
ducting gap and do not depend on the applied magnetic field. If we denote with E core 
the energy of the vortex core states measured with respect to the zero-bias voltage, then 



( 8 ) See Ref. [6] and references therein. 
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there is a simple relation between the state energy and the superconducting gap [6]: 

(4.33) E core ~ ± 0.3 A p , 

where the slope were obtained by a linear fit to all data points, and A p is the supercon- 
ducting gap. We have already noted that the data are consistent with the relation (see 
Fig. 8): 

(4.34) „ 4.3 , 

so that, using Eq. (4.2), we may rewrite Eq. (4.33) in terms of the pseudogap: 

(4.35) E core ~ ± 0.32 A 2 . 

To interpret this last remarkable relation, we refer to the Schrodinger equation for the 
d-wave hole bound states in presence of an external magnetic field, Eq. (4.18) with 
m = ±2. We have shown that, neglecting the Zeeman splitting, A±2(H, S) ~ A2(<5). 
Thus, we infer that the pseudogap without the external magnetic field fixes the zero-bias 
voltage in the core of an isolated vortex. Moreover, from Eq. (4.18) it follows that there 
are two states whose energies measured with respect to the pseudogap arc given by: 

(4.36) E core = ± , 

m h c 

where H core is the magnetic field inside the vortex core. To estimate H core , we note that 
the size of the isolated vortex can be determined by: 

(4-37) 

^A 2 

where ^ is the hole reduced mass. Thus, the magnetic field in the vortex core can be 
obtained from: 

(4.38) it & H core ~ (j)o , 

i.e. 

2irhc 1 
a 



(4.39) H core _ 

2e 7r£v 



Inserting Eqs. (4.39) and Eq. (4.37) into Eq. (4.36), we get: 
(4.40) E core ~ ± 1 A 2 



which is in satisfying agreement with Eq. (4.35). 
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4'4. Fermi arcs and quantum oscillations . - The normal-state properties in the un- 
dcrdopcd region of high-temperature cuprates are highly anomalous. It is now well- 
established by many angle- resolved photocmission studies ( 9 ) that low-energy excita- 
tions are characterized by Fermi arcs, namely truncate segments of a Fermi surface. 
The resistivity obeys p ~ T over a wide range of temperatures in striking contrast to a 
Fermi- liquid law p ~ T 2 , expected for a conventional metals. Moreover, several recent 
studies [9, 59, 60, 61, 62, 63] report unambiguous identification of quantum oscillations 
in high magnetic fields. Interestingly enough, the measured low oscillation frequencies 
reveals a Fermi surface made of small pockets. In fact, from the Luttinger's theorem 
and the Onsangcr relation [64, 65] between the frequency and the cross-sectional area 
of the orbit, it turns out that the area of the pocket correspond to about 2 — 3 % of 
the first Brillouin zone area in sharp contrast to that of overdoped cuprates where the 
frequency corresponds to a large hole Fermi surface [63]. In addition, the authors of 
Ref. [66] reported the observation of a negative Hall resistance in the magnetic-field 
-induced normal state, which reveals that these pockets are electron-like rather than 
hole-like. More recently [67], experiments on the second harmonic quantum oscillations 
in undcrdoped high-temperature cuprates lead to the conclusion that there exists only 
a single underlying quasi-two dimensional Fermi surface pocket. Moreover, it turns out 
that the pocket is most likely associated with states near the nodal region of the Brillouin 
zone. This nodal pocket comprises quasiclectron carrier to explain both the high field 
quantum oscillations with negative Hall and Seebeck [68] effect. Finally, the recent stud- 
ies of YBa2Cu3Oe.se at very high magnetic fields indicate that the specific heat exhibits 
the conventional temperature dependence and quantum oscillations expected for a Fermi 
liquid. Moreover, the magnetic field dependence of the quasiparticlc density of states 
follows the yH behavior pointing to a d-wave superconducting gap. Thus, the specific 
heat data demonstrate the surprising coexistence of the signature of a Fermi liquid and 
a d-wave superconducting gap over the entire magnetic field range measured. 
In the present Section we show that in our approach these features can be accounted 
for as a consequence of the pseudogap ^2(6). We have already seen that, in the un- 
dcrdoped region S < 5*, the holes are bound into pairs with k-space wavefunction given 
by Eq. (3.32). Even though the pair wavefunction vanishes along the nodal directions 
k x = ±k y (9k = ±f ), there are not nodal low-lying hole excitations. This is due to the 
fact that the pairing of the holes is in the real space and not in momentum space. On 
the other hand, due to the rotational symmetry, we may freely perform rigid rotations of 
the pairs without spending energy. The rigid rotations of pairs is equivalent to hopping 
of electrons according to the hopping term in the Hamiltonian Eq. (1.1): 



(4.41) 




<i,j>,o 



We may diagonalize this Hamiltonian by writing: 



(4.42) 




( 9 ) See Ref. [2] and references therein. 
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to get: 

(4.43) H ( Q e) = £ ef $t(k,<r)$ e (k,<r) 



(4.44) el e) = -2t [ cos (k x a ) + cos (k y a ) } 
In the small-k limit we have: 

(4.45) = Y,jS ti(k,<j)Mk,a), 



where: 

h 2 

(4.46) m* e = —3 ~ 2.17 m e . 

Since there are 1 — 8 electrons per Cu atoms, from the Hamiltonian Eq. (4.45) we may 
determine the electron Fermi energy: 



(4.47) 4 e) = . «o kf ~ VM^) 

2 m? 



At first glance, one expects that the quasielectrons fill in momentum space the circle 
with radius k^p) (the electron Fermi circle). However, one should keep in mind that the 
hopping of electrons is possible thanks to the paired holes. Since in momentum space the 
wavefunction of a given pair is spread over a region around k ~ ^ , we see that the wave- 
function of quasielectrons is likewise localized on a region in k-space around Thus, 
the quasielectrons do not have the needed coherence to propagate over large distances 
with a well defined momentum. However, this argument does not apply along the nodal 
directions where the momentum space hole-pair wavefunction vanishes. Therefore we 
are led to the conclusions that there are coherent quasielectrons that fill small circular 
sectors of the electron Fermi circle around the nodal directions k x = ±k y . Thanks to the 
rotational symmetry, we have four circular sectors with the same area. Since the number 
of coherent quasielectrons is determined by the doping fraction of holes (assuming that 
all the holes are paired), we obtain (see Fig. 11): 

(4.48) 4, = i2i \ { kf f 6fa , 
i.e. 

(4.49) Op a 



2 1-5 

We see, thus, that the Fermi surface is made by four Fermi arcs in agreement with the 
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Fig. 11. - The Fermi sector and Fermi arc in the first quadrant of the Brillouin zone. 



angle-resolved photoemission data. Moreover, the area of the Fermi sector with respect 
to the area of the Fermi circle in overdoped region turns out to be: 

(4.50) AfA - 1 6 



Aoverdoped 4 1 ~t~ & 

For the typical value of hole doping fraction S ~ 0.1, we infer from Eq. (4.50) that Afa 
is about 2.3 % of the first Brillouin zone area in the overdoped region, in satisfying agree- 
ment with quantum oscillation experimental results. Note that the effective quasielectron 
mass measured in quantum oscillation experiments lies in the range to* = (2.0 — 3.0) m e , 
in agreement with our estimate Eq. (4.46). 

The drastic reduction of the available states in momentum space leads to a peculiar 
quasielectron conductivity. To sec this, we may employ the simple Drude model [70] for 
the frequency-dependent conductivity: 



(4.51) a(u) ~ — ^— , a 



1 — iiUT TO* 

where, according to Eq. (4.48), n e ~ -4-. Since the nodal quasielectrons suffer almost 
unidimcnsional scatterings, we get for the relaxation time the quite general estimate: 



(4.52) j \ de' ~ e ~ k B T 
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and Eq. (4.51) becomes: 

1 1 k B T 

For fku <C fcsT we obtain the conductivity: 

n e e 2 h e 2 5 



(4.54) cr ~ cr 



On the other hand, if ftio 3> UbT we get: 

(To . n e e 2 . e 2 (5 1 



(4.55) o-(w) 



«o m e w 



These relations are in qualitative agreement with several observations. In particular 
Eq. (4.54) should be valid up to temperatures of the order of the pseudogap temperature 

Let us conclude this Section by discussing the peculiar dependence of the specific heat 
on the applied magnetic field reported in Ref. [69]. Firstly, we evaluate the Sommerfeld 
coefficient for H = 0. A standard calculation gives: 

(4,e, - ha ^ V J^ NAkl 

where Na is the Avogadro's number. The authors of Rcf. [69] found that at very high 
magnetic fields the specific heat in an underdoped high-temperature superconductors 
exhibits both the conventional temperature dependence and quantum oscillations ex- 
pected for a Fermi liquid. The oscillatory component of the specific heat is given by the 
Lifshitz-Koscvich formula [65]. On the other hand, for the Sommerfeld coefficient they 
report: 

(4.57) ls {H) = ls {H = 0) +A C ^H , 
with 

(4.58) j s (H = 0) = 1.85 ± 0.06 mJ mol' 1 ^ 2 , A c ~ 0.47 mJ mol^ 1 K~ 2 . 

To compare our Eq. (4.56) with Eq. (4.58) one should keep in mind that there are two 
CuO-2. planes per formula unit in YBCO. Moreover, the planar CuC>2 densities in high- 
temperature cuprates are difficult to estimate. Assuming 8 ~ 0.1 and using the effective 
electron mass Eq. (4.46), we find (in MKS units): 

(4.59) 7 S ~ 0.76 mJ mol~ 1 K~ 2 , 

which differs from measured value in Eq. (4.58) by about a factor of 2. As far as field- 
induced Sommerfeld coefficient is concerned, the behavior "f s (H) ~ \/~H is interpreted as 
evidence of the presence of a d-wave superconducting gap. In fact, it was pointed out [71] 
that the increase of the Sommerfeld coefficient with the external magnetic field is due 
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to the line nodes of the superconducting gap. The evidence of quantum oscillations 
and field-induced Sommerfeld coefficient in the electronic specific heat points to the 
surprising coexistence of a Fermi liquid and a d-wave superconductor. In our approach, 
we now show that the quasielectron nodal Fermi liquid leads naturally to the observed 
Sommerfeld coefficient j s (H) ~ \[H. In the mixed state, according to Ref. [71] the 
energies of the low-lying excitations circulating around a vortex are shifted by Dopplcr 
effect. Since it results that kp^o 1, we may employ the semiclassical approach which 
considers the momentum and position of quasiparticles as commuting variables. The 
effects of the supercurrent circulating around a vortex is accounted for by a Doppler 
shift of the quasiparticlc energy. Thus, for the quasielectrons we have: 



(4.60) 



s(k,r) ~ + hk-v s (r), 

2 m* 



where v s (r) is the supervelocity that, according to Eq. (4.15), is given by: 



(4.61) 



v s {r) 



h 6 r 



4 m* 



The Doppler shift Eq. (4.60) has a sizable effect on the density of state. To see this, wc 
need to evaluate the density of state at the Fermi surface averaged over the vortex: 



(4.62) 

where 



Af(0) 



dk 



dr 5 



e(k, r) ~ e 



1 



irR 2 ' 



(4.63) 



is the average distance between vortices. We obtain: 
(4.64) 



dOi 



FA 



dk k / dO. 



dr 



So 



irR 2 



h 2 k 2 h 2 k sin (9 r 



2m* 



Am*,r 



y k) _ Je) 



where the angular integration on 6^ is performed over the four Fermi arcs (FA). The 
integration over k can be done easily. One obtains: 



(4.65) Af(0) ~ -L 



47r 2 h 2 m* h j FA 



dOi. / dO. 



R dr 



ttR 2 



Finally the angular integrations give: 



(4.66) 
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or, using Eqs. (4.63) and (4.49): 

V2 S {m* e 



(4.67) A/"(0) 



4^ 1-5 K* m * h k f R 




Since the high-temperature superconductors are extreme type-II superconductors, one 
should average the density of state over the vortex distribution. However, from Eq. (4.67) 
it is evident that the quasiparticle energy Doppler shift induced by vortices leads to the 
expected field-induced Sommerfeld coefficient ^ S {H) ~ \fH. 

5. The physics of the over and optimal doped regions 

In the overdoped region we have already stressed that the normal state properties 
can be described reasonably well by the Fermi liquid picture, although still with some 
electronic correlations, with a hole density corresponding to 1 + 5 holes per Cu atom. 
Indeed, as long as S > 5*, the effective two-body potential Eq. (2.7) can be dealt with 
as a small perturbation. Therefore, the dynamics of the copper-oxide layers is governed 
essentially by our effective Hamiltonian Hq, Eq. (2.18). Note that the correlation effects 
due to the on-site repulsion U are built into Hq by means of the effective hole mass m* h . 
It is interesting to note that: 

mi U 

(5.1) — ^ = — ~ 2.5 . 

y ' m* At 

In fact, several quantum oscillation experiments give consistently an effective quasipar- 
ticle mass of carrier in the overdoped region which is approximately twice the effective 
quasiparticle mass of carriers in the underdoped region. 

It is interesting to compare the low-temperature linear term in the electronic specific heat 
of the quasiholes with the one of the nodal quasielectrons evaluated in Sect. 4'4. Since 
the quasiholes behave like a two-dimensional normal Fermi liquid, we obtain instead of 
Eq. (4.56): 

(5.2) 7s * 3 N A k% . 

In fact, it is easy to check that the Sommerfeld coefficient given by Eq. (5.2) exceeds by 
about an order of magnitude the quasielectron Sommerfeld coefficient. Using our hole 
effective mass we obtain (MKS units): 

(5.3) j s ~ 8.53 mJ mol~ x K~ 2 . 

Our determination of 7 S compares quite well with the observed value: 

(5.4) 7s ~ 7 ± lmJmor'r 2 

as obtained from direct measurements in polycrystalline Tl2Ba2CuOe+s [72], indepen- 
dently on 6, and also from indirect determination through the effective hole mass = 
(5.2±0.4)m e in Ref. [63]. 
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Concerning the superconducting phase in the overdoped region, we already argued that 
this can be described within the conventional BCS framework. In fact, there is a large 
body of evidence that the weak-coupling d-wavc BCS approach accounts for many of the 
low-energy and low-temperature properties of overdoped copper oxides. Here, we restrict 
ourself to the discussion of a clear signature of d-wavc BCS pairing in momentum space, 
namely the ratio of the BCS gap to the critical temperature anticipated in Eq. (3.61): 



(5.5) 



A 



BCS 



2 7T 



k B T c 



2.14 



To derive Eq. (5.5), we observe that the gap equation Eq. (3.47) at finite temperature 
reads: 



/ \ »,a 1 V(k-k')A(k') 

(5.6) A(k) = -- / : tanh 

£| + |A(fc')l 2 



q, + |A(fc 



AI2 



2k B T 



Thus, instead of Eq. (3.56) we have: 
(5.7) 



1 



-Vn 



dk' 



[cos {20 k ,)f 



tanh 



[A 



BCS 



cos(20 fe ,)] 5 



[A BCS cos(2e fc ,)F 



2k B T 



Since the critical temperature is defined by A B cs(T c ) = 0, we obtain: 



(5. 



K 2 



JL 

(2ny 



2n in [cos(20)] 2 

dO v n tanh 

Isl 



2k B T n 



i.e. 



(5.9) 



V2 m* h 
2n h 2 



f ^ tanh 




10 S 


2k B T c 



Upon using Eq. (3.60) we may rewrite Eq. (5.9) as 
(5.10) 1 



f2k B T c dz , . . 

A 2 / — tanh [z\ 

Jo z 



A 2 



which leads to: 
(5.11) 



k B T c ~ y exp(-ln^+7) exp[-^-] 



Combining Eq. (5.11) with Eq. (3.59) one obtains Eq. (5.5). In Figure 12 we display the 
doping dependence of the ratio Eq. (5.5) using the data reported in Fig. 4 of Ref. [73]. 
It is evident that in the overdoped region 6 > 0.20 the gap values approach closely onto 
the theoretical prediction of the weak-coupling d-wave BCS superconductivity. 
Let us now discuss the optimal doped region 5 ~ 5 op t. According to our picture, the 
low-doping region is characterized by a large pscudogap A 2 ((5) which decreases rapidly 
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H 

< 




0.35 



Fig. 12. - Doping dependence of the ratio of superconducting gap A over the critical temperature 
T c . The data has been extracted from Fig. 4 of Ref. [73]. The dashed line is the weak coupling 
d-wave BCS value rAp- ~ 2.14. 



with increasing doping, while the superconductivity is driven by the finite-temperature 
Berezinskii-Kosterlitz-Thouless transition with a critical temperature that increases with 
doping. On the other hand, in the overdoped region the superconducting transition is 
well described by the conventional weak-coupling d-wave BCS pairing in momentum 
space. The decrease of the critical temperature with increasing doping level originates 
from the underlying reduction in the pairing strength due to the disturbing of the anti- 
ferromagnetic correlations with overdoping. From our Fig. 7 we see that the two critical 
temperatures meet in the optimal doped region 5 ~ 5 op t, explaining in a natural way 
why the critical superconducting temperature reaches here the maximum value. Since 
the Berezinskii-Kosterlitz-Thouless temperature is comparable to the BCS critical tem- 
perature, this makes fluctuations much more important. Even more, Fig. 7 suggests that 
also the pseudogap temperature T* becomes comparable to both Tb-k-t an d Tbcs- 
We conclude, thus, that in this region of the phase diagram there are competing phases 
and it is quite difficult to precisely pin down the underlying physics. On the other hand, 
as our Fig. 7 is suggesting, there is a large body of evidence that pseudogap energy scale 
does not extend into the heavily overdoped region of the phase diagram, but rather col- 
lapse at a well defined critical concentration around 8* ~ 0.20. We would like to stress 
that the presence of an expanded phase-fluctuation regime in the optimal doped region 
may explain why certain spectroscopic probes, like angle-resolved photocmissions, tend 
to advocate scenarios in which the pseudogap temperature tracks the superconducting 
dome rather than vanishing inside it. The scenario advanced in this paper is supported by 
the results presented in Ref. [74] . By means of a systematic electronic Raman scattering 
study of a mercury-based single layer cuprate as a function of doping level, the authors of 
Ref. [74] revealed the existence of a breakpoint close to optimal eloping, below which the 
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antinodal gap (the pseudogap) is gradually disconnected from superconductivity. The 
nature of both the superconducting and normal state is distinctly different on each site 
of this breakpoint. This is consistent also with electronic specific heat measurements. 
In fact, it turns out that [75, 76] in optimally and overdoped samples the Sommerfeld 
coefficient is almost constant for T > T c , in accord with our Eq. (5.2), while for under- 
doped samples there is considerable drop in the Sommerfeld coefficient for T < T* (see 
our Eq. (4.56)). Moreover, there is evidence for carrier density which decreases smoothly 
from in the overdoped phase to -4- in the underdopcd phase. In fact, in Rcfs. [77, 781 
it is shown that it exists a critical point 8* where the pseudogap temperature T* goes 
to zero and which lies in the superconducting dome. Just below 8* there is a crossover 
from a metal with large hole-like Fermi surface, consistent with rih ~ ^r, to a metal 

a 

with low density of charge carrier, which from Hall and Seebeck effect measurements are 

quasiclcctrons, with n e ~ -4-. 

a ° 

Before concluding this Section, wc would like to remark that the optimal doped region 
of hole-doped cuprates is the most difficult to precisely characterize. In fact, according 
to our current theoretical understanding, in this region there is competition between the 
BCS gap and the pseudogap, and the presence of sizable phase fluctuations strongly limit 
our mean-field treatment. Nevertheless, it is worthwhile to disciss the dependentee of 
the maximum critical temperature T™ ax on the microscopic parameters of the model. 
As already discussed we have: 

(5.12) r c nax ^ T BCS (6 opt ~ 8*) . 

Now, from Eqs. (2.5), (2.9), (2.17), (3.46) and (3.60) we obtain: 

where we recall that 5* is the hole doping fraction where the pseudogap A 2 vanishes. 

Apparently Eq. (5.13) suggests that one can increase T™ ax by increasing jj- and/or 

.2 

decreasing 8*. However, it turns out that the increase of jj leads also to an increase of 
8* and vice versa, such that the two effects compensate in Eq. (5.13). This explains why 
the optimal hole doping fraction is very close to 8* and attains the universal value: 

(5.14) 8 opt ~ 8* ~ 0.20 

in the whole family of cuprates. 

6. The electron doped cuprates 

The vast majority of high-temperature superconductors are hole-doped. Nevertheless, 
in this Section we will touch upon high-temperature electron-doped cuprates (see Ref. [19] 
for a review). The most evident difference between electron and hole doped cuprates is 
that the antiferromagnetic phase at 8 = is much more robust in the electron-doped 
material and it persists to much higher doping levels. If we adopt the effective single- 
band Hubbard model of Sect. 1 for the microscopic description of the CuOi layers, then 
we see that the electrons injected into the copper-oxide plane would behave like the 
holes in the overdoped region of the hole-doped cuprates. However, due to the stronger 
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antiferromagnetic correlations in the ground state for 5 > 0, the electron carrier density 
is: 

(6.1) n e ~ \ , 



since one electron of the (1 + 6) electrons per Cu atom is frozen in the half- filled Hubbard 
band antiferromagnetic configuration. In our approach the effective Hamiltonian for the 
injected electrons should be given by Eqs. (2.18) and (2.6). In particular, the electron 
effective mass agrees with Eq. (2.17): 



(6.2) ml ~ — — ~ 5.41 m e 

with the previously stated values of the model parameters. However, due to the increased 
role played by the fermion correlations in the antiferromagnetic background, we expect 
that the effective two-body potential in the interaction Hamiltonian Eq. (2.6) is reduced 
to some extent. To take into account this effect, we shall assume that the effective 
potential V(f\ — ra) has the same form as in Eq. (2.7) but with a somewhat reduced 
interaction length: 



(6.3) r ( e) (6) = 4.15 a Q f 1 - — J , 6 C ^ 0.35. 

We now show that this small variation in the range of the effective two-body interaction 
potential leads to a dramatic change in the phase diagram of the electron-doped cupratcs. 
In fact, by solving the real-space d-wave bound state Schrodinger equation Eq. (3.11) we 
find a drastically reduced pseudogap (sec Fig. 13): 

A (e) f<5 = 0) 

(6.4) T? e) (5 = 0) = ~ 46 X , 6* ~ 0.05 . 

y ' 2 Kb 

This means that, according to the discussion in Sect. 3'2, the pseudogap superconduct- 
ing phase of the undcrdoped region of hole-doped cupratcs is absent. Therefore, the 
superconducting transition in electron-doped cuprates is driven by the weak-coupling 
d-wavc BCS pairing in momentum space. In this case the critical BCS temperature is 
given by Eq. (3.62): 

(6.5) k B T B % s (5) ~ — A< e) (0) exp[ 



with A^ e) (0) given by Eq. (6.4), 



(6.6) A< e) = ^ , 



„0)„ 
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Fig. 13. - The pseudogap and BCS critical temperatures (in Kelvin) versus the doping fraction 
for the electron-doped cuprates. 



and 



(6.7) 



„(«0 



ao 



The BCS critical temperature Eq. (6.5) as a function of the electron doping level is 
displayed in Fig. 13. From the schematic phase diagram we infer that the maximal critical 
temperature is strongly reduced with respect to the hole-doped cuprates. Moreover, the 
superconducting phase manifests itself in a rather narrow interval of doping around S ~ 
0.12. These general aspects of our phase diagram seem to be in qualitative agreement with 
observations. However, a thorough understanding of the superconductivity in electron 
doped cuprates must await more detailed experimental informations. 



7. Summary and Conclusions 

The driving principle of the approach presented in this paper has been that the high- 
temperature superconductivity can be understood within some framework along the line 
of the microscopic theory of Bardcen, Cooper, and Schricffer. The results presented in 
this paper summarize our attempts over the last years to construct a comprehensive 
model that gives all aspects of the unusual behavior seen in the various regions of the 
phase diagram of the cuprates. It is obvious that we must rely heavily on some assump- 
tions. First, we assumed that the physics of the high-temperature cuprates is deeply 
rooted in the copper-oxide planes. This allowed us to completely neglect the motion 
along the direction perpendicular to the CuC>2 planes. Our second main assumption has 
been that the single-band effective Hubbard model is sufficient to account for all the 
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essential physics of the copper-oxide planes. Both assumptions are well accepted and 
superbly illustrated in the Anderson's book [11]. Accordingly, we have proposed an ef- 
fective Hamiltonian that is able to govern the dynamics of the electrons or holes injected 
into the undoped copper-oxide planes. We arrived at our effective Hamiltonian by using 
arguments well-known in the literature on the motion of charge carriers in an antifcr- 
romagnetic background. Notwithstanding, we were unable to offer a truly microscopic 
derivation of the effective Hamiltonian. So that our arguments to determine the effective 
Hamiltonian, albeit suggestive, cannot be considered a first principle derivation. In spite 
of that, we showed that the effective Hamiltonian offers us a consistent picture of the 
high transition temperature cuprate superconductors. Due to the reduced dimensionality 
the two-body attractive potential admits real-space d-wave bound states. The binding 
energy of these bound states, which plays the role of the pscudogap, decreases with in- 
creasing doping until it vanishes at a certain critical doping 6 = 5*. This allows us to 
reach the conclusion that the key features of the underdoped side of the phase diagram 
arc controlled by very strong pairing that is phase-disordered by thermal fluctuations. 
We also were able to estimate the upper critical magnetic field whose strength and doping 
dependence are in good agreement with observations. We obtain a clear explanation for 
the so-called vortex region and for the puzzling discrete states bounded to the core of 
isolated vortices. We have discussed how the presence of the pseudogap is responsible 
for the formation of the quasielectron nodal Fermi liquid which, in turn, leads to the 
Fermi arcs observed in angle-resolved photoemission studies and to the Fermi pockets in 
quantum oscillation experiments. We have found a natural explanation for the peculiar 
dependence of the specific heat on the applied magnetic field observed in underdoped 
cupratcs. However, we would like to stress that we did not attempt a complete discussion 
of the strange behavior of the cuprates in this region. In particular, as already pointed 
out in Sect. 4, we did not discuss the so-called nodal gap. In fact, we plan to present in 
a separate paper a detailed discussion on the origin and temperature dependence of the 
nodal gap, on the thermal conductivity, on the departure from the Wiedemann-Franz 
law and on the London penetration depth. 

In our model the ovcrdoped region is realized for hole doping in excess of the critical 
doping 5* where the pseudogap vanishes. In this region the conventional d-wave BCS 
framework account for many of the low-energy and low-temperature properties of the 
copper oxides, in agreement with several observations. Finally, we explain naturally why 
the superconducting critical temperature reaches its maximum in the optimal doped re- 
gion. We pointed out that in this region the competition between the pseudogap and 
the d-wave BCS gap together with the enhanced role of the phase fluctuations makes 
questionable the usually adopted mean-field approximation. We have also briefly dis- 
cussed the comparatively less-studied electron-doped cuprates. We suggest that the 
dome-shaped superconducting region in electron-doped copper-oxides can be described 
by the conventional weak-coupling d-wave pairing in momentum space. 

In conclusion, the fact that a relatively simple model for the effective Hamiltonian 
allowed to recover several observational features of the high-temperature copper oxides, 
suggests that the approach presented in this paper could be considered a significative step 
toward a microscopic explanation of the high-temperature superconductivity in cuprates. 

REFERENCES 

[1] J. G. Bednorz and K. A. Miiller, Z. Phys. B 64 (1986) 189 



THE HIGH TEMPERATURE SUPERCONDUCTIVITY IN CUPRATES 



39 



[2] A. Damascelli, Z. Hussain and Z. X. Shen, Rev. Mod. Phys. 75 (2003) 473 

[3] G. Deutscher, Rev. Mod. Phys. 77 (2005) 109 

[4] D. N. Besov and T. Timusk, , Rev. Mod. Phys. 77 (2005) 721 

[5] P. A. Lee, N. Nagaosa, and X. -G. Wen, Rev. Mod. Phys. 78 (2006) 17 

[6] 0. Fischer, H. Kugler, I. Maggio-Aprile and C. Berthod, Rev. Mod. Phys. 79 (2007) 353 

[7] P. A. Lee , Rep. Prog. Phys. 71 (2008) 012501 

[8] S. Hiifner, M. A. Hossain, A. Damascelli and G. A. Sarvatzky, Rep. Prog. Phys. 71 (2008) 
062501 

[9] S. E. Sebastian, N. Harrison and G. G. Lonzarich, Rep. Prog. Phys. 75 (2012) 102501 
[10] P. W. Anderson, Science 235 (1987) 1196 

[11] P. W. Anderson, The Theory of Superconductivity in the High-T c Cuprates, Princeton 

University Press, Princeton, New Jersey (1997) 
[12] P. W. Anderson, Phys. Rev. 115 (1959) 2 
[13] D. A. Bonn, Nature Physics 2 (2006) 159 
[14] Ch. Niedermayer et al, Phys. Rev. Lett. 80 (1998) 3843 
[15] S. -W. Cheong et al, Phys. Rev. Lett. 67 (1991) 1791 

[16] T. E. Mason, G. Aeppli and H. A. Mook, Phys. Rev. Lett. 68 (1982) 1414 
[17] A. N. Lavrov, L. P. Kozeeva, M. R. Trunin and V. N. Zverev, Phys. Rev. B 79 (2009) 
214523 

[18] J. Bardeen, L. N. Cooper and J. R. Schrieffer, Phys. Rev. 106 (1957) 162; Phys. Rev. 108 
(1957) 1175 

[19] N. P. Armitage, P. Fournier and R. L. Greene, Rev. Mod. Phys. 82 (2010) 2421 

[20] A. P. Balachandran, E. Ercolessi, and G. Morandi, Hubbard Model and Anyon 

Superconductivity, World Scientific Publishing Co. Inc. (1990) 
[21] K. Huang and E. Manousakis, Phys. Rev. B 36 (1987) 8302 
[22] J. E. Hirsch, Phys. Rev. Lett. 59 (1987) 228 
[23] S. A. Trugman, Phys. Rev. B 37 (1988) 1597 

[24] E. W. Carlson, V. J. Emery, S. A. Kivelson and D. Orgad, Concepts in High Temperature 

Superconductivity , cond-mat/0206217 
[25] J. R. Schrieffer, X. -G. Wen, and S. -C. Zhang, Phys. Rev. Lett. 60 (1988) 944; Phys. Rev. 

B 39 (1989) 11663 
[26] L. N. Cooper, Phys. Rev. 104 (1956) 1189 
[27] M. R. Schafroth, Phys. Rev. 100 (1955) 463 
[28] J. M. Blatt and S. T. Butler, Phys. Rev. 100 (1955) 476 
[29] P. W. Anderson, Rev. Mod. Phys. 38 (1966) 298 

[30] A. J. Legget, Quantum Liquids, Oxford University Press, Oxford, UK (2008) 

[31] V. L. Berezinskii, Sov. Phys. JETP 32 (1971) 493; Sov. Phys. JETP 34 (1972) 610 

[32] J. M. Kosterlitz and D. J. Thouless, J. Phys. C 6 (1973) 1181 

[33] B. Vignolle et ai, Nature 455 (2008) 952 

[34] P. W. Anderson and P. Morel, Phys. Rev. 123 (1961) 1911 

[35] H. Won and K. Maki, Phys. Rev. B 49 (1994) 1397 

[36] T. Timusk and B. Statt, Rep. Prog. Phys. 62 (1999) 61 

[37] J. L. Tallon and J. W. Loran, Physica C 349 (2001) 53 

[38] U. Chatterjee et ai, Nature Physics 6 (2010) 99 

[39] K. Tanaka et al, Science 314 (2006) 1910 

[40] R.-H. He et al, Nature Physics 5 (2009) 119 

[41] A. Pushp et al, Science 324 (2009) 1689 

[42] M. Hashimoto et al, Nature Physics 6 (2010) 414 

[43] T. J. Reber et al, Nature Physics 8 (2012) 606 

[44] S. Sakai et al, Exploring the Dark Side of Cuprate Superconductors: s-wave Symmetry of 

the Pseudogap, arXiv: 1207.5070 
[45] T. Yoshida et al, Coexisting pseudo-gap and superconducting gap in the high-T c 

superconductor La2- x Sr x CuOi, arXiv: 1208.2903 
[46] M. Hashimoto et al, Phys. Rev. B 86 (2012) 094504 
[47] M. Kugler et al, Phys. Rev. Lett. 86 (2001) 4911 



40 



PAOLO CEA 



[48] I. Hetel, T. R. Lemberger and M. Randeira, Nature Physics 3 (2007) 700 

[49] J. Corson et al, Nature 398 (1999) 221 

[50] V. J. Emery and S. A. Kilvelson, Nature 374 (1985) 434 

[51] M. Tinkham, Introduction to Superconductivity, Second Edition, Mc Graw-Hill, Inc., New 
York (1996) 

[52] Y. Wang, L. Li and N. R Ong, Phys. Rev. B 73 (2006) 024510 

[53] L. Li et al, Nature Physics 3 (2007) 311 

[54] J. Chang et al, Nature Physics 8 (2012) 751 

[55] Z. A. Xu et al, Nature 406 (2000) 486 

[56] C. Capan et al, Phys. Rev. Lett. 88 (2002) 056601 

[57] L. Li, Y. Wang, J. G. Checkelsky and M. J. Naughton, Physica C 460-462 (2007) 48 

[58] C. Caroli, P. G. de Gennes and J. Matricon, Phys. Lett. 9 (1964) 307 

[59] N. Doiron-Leyrand et al, Nature 447 (2007) 565 

[60] E. A. Yelland et al, Phys. Rev. Lett. 100 (2008) 047003 

[61] A. F. Bangura et al, Phys. Rev. Lett. 100 (2008) 047004 

[62] C. Jaudet et al, Physica B 404 (2009) 354 

[63] A. F. Bangura et al, Phys. Rev. B 82 (2010) 140501 

[64] A. A. Abrikosov, Introduction to Theory of Normal Metal, Academic Press, New York 
and London (1972) 

[65] D. Shoenberg, Magnetic Oscillations in Metals, Cambridge University Press, Cambridge 
(1984) 

[66] D. Le Boeuf et al, Nature 450 (2007) 533 

[67] S. Sebastian et al, Nature Communications 2 (2011) 471 

[68] F. Laliberte et al, Nature Communications 2 (2011) 432 

[69] S. C. Riggs et al, Nature Physics 7 (2011) 332 

[70] N. W. Ashcroft and N. D. Mermin, Solid State Physics, Harcourt College Publishers (1976) 

[71] G. E. Volovik, Sov. Phys. JETP 58 (1993) 469 

[72] J. M. Wade et al, J. Supercond. 7 (1994) 261 

[73] Y. Wang et al, Phys. Rev. B 76 (2007) 064512 

[74] W. Guyard et al, Phys. Rev. B 77 (2008) 024524 

[75] J. W. Loram, K. A. Mirza, J. R. Cooper and J. L. Tallon, Physica C 282-287 (1997) 1405 

[76] J. W. Loram, J. Luo, J. R. Cooper and J. L. Tallon, J. Phys. Chem. Sol. 62 (2001) 59 

[77] R. Daou et al, Nature Physics 5 (2009) 31 

[78] R. Daou et al, Phys. Rev. B 79 (2009) 180505 



